A new paradigm for two-dimensional spin-liquids 
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Motivated by the geometry of the materials Na2Ti2As20 and Na2Ti2Sb20, we study a square-lattice 
Heisenberg antiferromagnet, with spins located at the bond-centers. The largest exchange constant 
J couples neighboring spins in a given row or column. This leads to a mesh of isolated spin-chains 
running along the X and Y axes. A weaker exchange constant J' couples the nearest-neighbor spins 
on the lattice. Classically, J' fails to fix the relative spin orientation for different chains and hence the 
ground state is highly degenerate. Quantum order by disorder effect is studied by spin-wave theory 
and numerical methods. It is shown that a 4-sublattice order is favored by quantum fluctuations. 
However, several arguments are presented that suggest that the ground state of the system remains 
disordered, thus providing us with a paradigm for a two-dimensional spin-liquid. 

PACS numbers: 75.I0.Jm, 75.40. Cx, 75.40. Gb, 75.50.Ee 



In recent years many new materials have been dis- 
covered which exhibit novel magnetic behavior. Var- 
ious aspects of quantum magnetism, including quan- 
tum critical phenomena and existence of spin-disordered 
ground states with spin-gaps have been observed in a 
variety of Cuprates, Germanates, Vanadates and other 
low-dimensional materials. One interesting fact that has 
come to light in these studies is that the geometrical 
arrangement of the transition metal and Oxygen ions 
can have a dramatic impact on the underlying micro- 
scopic spin-Hamiltonian and hence on the macroscopic 
magnetic properties of the system. For example in vari- 
ous Cuprates, the Cu-O-Cu bond angle is crucial for de- 
termining the effective exchange constant between Cop- 
per spins. Thus, Strontium Cuprates with certain sto- 
ichiometry behave as virtually decoupled spin-ladders 

, even though the separation of the Copper spins be- 
tween neighboring ladders maybe smaller than their sep- 
aration within a given ladder. In the CuGe03 B and 
CaV n 2 „ + i it is also believed that superexchangc 
between the transition metal ions is mediated by Oxy- 
gen and because of the geometry of various bond angles 
and occupied orbitals the second neighbor interactions 
are substantial compared to nearest neighbor ones. This 
leads to various interesting quantum phase transitions 
and spin-gap behavior in these materials. 

Here, we consider a Heisenberg model with spins at the 
bond-centers of a square lattice: 



H — J ^ Si ■ Sj + J 1 S Sj ■ Sj, 



(1) 



<«,j> 



<i,3> 



with J' << J. The interactions are shown in Fig. 1. 
The exchange J couples neighboring spins in a given 
row or column, whereas J' is the nearest neighbor cou- 
pling between rows and columns. In the absence of J' 
the system consists of a square mesh of decoupled spin- 
chains running along the X and Y axes. The motiva- 
tion for studying such a model comes from the materials 



Na 2 Ti 2 Sb 2 ( and also Na 2 Ti 2 As20) These lay- 

ered Titanates consist of planes of (Ti 2 Sb 2 0) 2 ~, where 
the Oxygen atoms form a square-lattice and the Tita- 
nium atoms sit at the bond centers of the lattice. The 
Antimony atoms sit above and below the centers of the 
elementary squares. The Titanium atoms carry spin-half. 
It is evident from the geometry that if the dominant ex- 
change interaction is provided by a direct overlap be- 
tween the Titanium orbitals, this material would behave 
as a nearest-neighbor Heisenberg model. If on the other 
hand the dominant interaction is mediated by Oxygen 
orbitals it would result in the largest interaction between 
neighboring spins in a given row or column, thus lead- 
ing to our Eq. [I] Q. Finally, if the dominant exchange 
is mediated by the p-orbitals of Antimony, the Hamilto- 
nian in Eq. [l] could still be appropriate, the lattice being 
dual to the one shown. The magnetic properties of this 
material are currently under investigation ||. The rest 
of the paper will be devoted to a study of this model 
Hamiltonian. 
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FIG. f . The exchange interactions J and J' in our model, 
shown by solid lines and broken lines respectively. 
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We begin by studying this model in the linear spin- 
wave approximation. As finding the correct local spin 
configuration is based on energetic considerations, we ex- 
pect the spin-wave theory to be at least qualitatively cor- 
rect. Notice that the problem at hand has 'double' degen- 
eracy in the classical limit: it is evidently degenerate with 
respect to relative angle 9 between quantization axes on 
vertical and horizontal chains, as well as with respect to 
relative orientation 6 of quantization axes on the neigh- 
boring parallel chains. Based on the fact that Hamilto- 
nian (|lj) classically has no spiral ordering, and on insight 
from a somewhat analogous problem of coupled planes 
, we restrict variation in 6 to two possible values: 
6 = (ferromagnetic ordering of neighboring parallel 
chains) and 6 — ir (antiferromagnetic one). This degree 
of freedom is then naturally represented by Ising-type 
discrete variable t^+i = 1(6 = 0) or r^+i = 0(6 = n), 
defined for each pair of chains (i, i + 1). The calculations 
are simplified greatly by choosing quantization axes on 
all sites such that the ordering is ferromagnetic. This is 
achieved by a unitary transformation parametrized by an 
angle 9 and a set of (r^+i). 

We now integrate out spins on the horizontal chains 
to find an effective Hamiltonian for the remaining spins. 
This is achieved by representing spins on vertical (hori- 
zontal) chains in terms of Holstein-Primakoff bosons a (6), 
and doing perturbation expansion in J'/ J. First nonzero 
contribution is of the type J' 2 e' S" +e (a)Sj +e ,(a) < 

Sf (b)Sj(b) >b, where averaging is with respect to the 
bare Hamiltonian of horizontal chains and <Si+ e spins be- 
long to vertical chains, i.e. J' acts along the < i, i + e > 
links. The average is nonzero when i and j belong to 
the same horizontal chain and a = f3. It is seen that 
effective interaction between spins on vertical chains is 
time-dependent, but strongest contribution comes from 
static (w n = 0) part of the average. Denoting J^_(i—j) = 
J' 2 < Sf(b)S f j(b) > b and J^(r) = J^(r - 1) + JJ(r) we 
find that geometry dictates the following coupling be- 
tween remaining parallel chains n and m, 

(2) 

Notice that long-range-ordered part of < Sf(b)S^(b) 
if any, cancels out due to AFM correlations along the 
chains, and the effective interaction decays at least as 
(n — m)~ 2 . It is also important to realize that J" is 
anisotropic in spin space, in particular J^ z ^ ~ cos 2 9. 
Notice that classically the system remains degenerate be- 
cause < 25^ i + S^ i+1 + S^^-l >= 0. Given the fast 
decay of the induced inter-chain interaction with the dis- 
tance we restrict ourselves to the strongest interaction, 
J"(l), between neareast chains. To find corrections to 
ground state energy due to this coupling we next treat 



J± as a perturbation (see Ref. [|T(| for description of the 
procedure) . After long algebra one finds quantum correc- 
tion to the ground state energy (per chain) of M parallel 
chains of length TV: 

M 

SE G .s. =-j-Y1 <%W/Jf cos 2 0^(1- 2r M+1 ), 

q v i=l 

(3) 

where e(q) = 2JS\ smq\ is the single AFM chain disper- 
sion in spin- wave theory, C = (< SqSq > + < SfiSf > 
) 2 /64, and explicit form of Jj_was used. Note that correc- 
tion is down by l/S factor, showing its quantum origin. 
Thus the energy is minimized by choosing all Ti^+i = 
(AFM configuration) and cos 2 9 = 1. This result agrees 
with well-known tendency of quantum fluctuations to fa- 
vor colinear structures ■ 

Having found the ground state configuration of spins, 
we can determine the spin-wave dispersion and calculate 
the reduction in sublattice magnetization due to quan- 
tum fluctuations. This reduction turns out to be diver- 
gent due to a remaining artificial degeneracy of the spin- 
wave spectra. Similar divergency arises in a quantum 
order by disorder calculations on other systems and is 
known to be removed by the higher-order quantum cor- 
rections to the spin- wave spectra ]T^] . This technically 
difficult calculation has not been done. We would like 
to stress that the ground state structure found in (|^) 
is determined by short-range spin correlations (i.e. by 
the correlations within the correlation length range) for 
which spin-wave approximation should work even if the 
sublattice magnetization vanishes. 

Let us now take a closer look at the elementary exci- 
tations of the single chain. It is convenient to perform 
"staggering" of the spin configuration so that Neel or- 
dering along the chain corresponds to the ferromagnetic 
one in the new representation. In this representation 
an elementary excitation of the chain is a domain wall 
(spinon) separating ferromagnetic domains of different 
orientation. The energy of a single domain wall is J z . 
Consider now unfrustrated coupling J± between nearest 
spins on neighboring parallel chains: however small Jj_ 
is, it immediately leads to the suppression of spinons be- 
cause the energy costs is proportional to J± x (length of 
the domain) and diverges in the thermodynamic limit. 
This is an intuitive reason for the stabilization of LRO 
in the system of unfrustrated coupled chains [||. But 
this is not true for our Hamiltonian (||), where each spin 
is coupled to the zero-spin combination of spins on the 
neighboring chains. Thus domain wall excitations of the 
single chain seems to survive in the presence of nonzero 
J', hinting to the possibility of the absence of LRO. 

The effect of the "mixing" interchain coupling J' can 
also be taken into account in the disordered phase, 
where the interspin interaction is isotropic. We write 
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the partition function of the Hamiltonian (1) in the in- 
teraction representation, where independent chains are 
considered as unperturbed system and interchain cou- 
pling J' as a perturbation fll4| . Performing trace over 
the horizontal chains first one finds effective interac- 
tion between spins on the vertical chains proportional to 
(J / (k)) 2 G^{k x ,uj n )S a (k,uj n )S a (k,uj n ) 1 similar to our 
previous spin-wave calculations. Here G^ (k x ,uj n ) is 
the rotationally-invariant one-dimensional spin Green's 
function in Matsubara representation, and J'(k) = 



J' cos % cos • 



Decoupling spins on different chains 



via a Hubbard- Stratonovich transformation |14[], the 2D 
Green's function of the spins on vertical chains becomes, 



G(k x , k y , u> n ) 



G^(k y ,UJ n ) 



l-(J'(k))2G(o)(fc ;E , W „)G(o)(fc !y ,^ n )' 



(4) 



The known form of G^(k x ,Lu n ) |l5| implies that (i) 
the singularity at ko = (n, tt) remains unchanged be- 
cause J'(ko) = 0, (ii) the uniform susceptibility is 
X = Xo/[l - (XoJ') 2 }, where Xo = VO 1 " 2 - 7 ) is suscepti- 
bility of single chain, and (iii) interchain coupling comes 
into play below T* ~ J' . Within this approximation, the 
two-dimensional system remains in the disordered criti- 
cal state at T = 0, and in the quantum-critical one at 
finite temperatures. This result should also be viewed 
with some caution due to the large degeneracies in the 
underlying problem which we do not know yet how to 
take into account completely. 

We now turn to numerical studies. Because of multi- 
ple energy scales, the development of short-range order at 
high temperatures is not necessarily indicative of the or- 
der for T ~ J 1 < < J. Thus we need to study the ordering 
tendency directly at low temperatures. As J 1 / J « 1, 
the coupling between the local order parameters on the 
neighboring parallel chains can be gotten by perturba- 
tion theory. We study the following correlation function 
between neighboring vertical chains: 



Cl2 



=<£(-i)<s M .£(-iys 2ii >, 



(5) 



here Sx,i represents spins on first vertical chain and S2,i 
represents spins on second vertical chain. A value of 
C12 > would imply ferromagnetic ordering whereas 
C12 < would imply antifcrromagnetic ordering. The 
leading order contribution to c\i requires the interactions 
to be mediated by at least two horizontal chains. To make 
the perturbation theory numerically tractable, we confine 
ourselves to finite chains along vertical and horizontal di- 
rections ( with periodic boundary conditions). For four 
spins in each of the chains, we can carry out perturbation 
theory using series expansion methods p2[ | . We find that 
C12 f» —9. x 10~ 6 (J'/J) 6 , which is negative confirming 
the tendency of the parallel chains to align antiferromag- 



netically. We have also used Lanczos methods to calcu- 
late this correlation function with upto 6-spins in each 
chain for \J'/J\ < 0.3. In all cases the coupling remains 
antiferromagnetic and very weak. For the 6-spin system 
C12 appears to scale as (J'/ J) 4 as expected from general 
arguments and spin-wave theory presented earlier. Com- 
plete diagonalization of a 32-site system is currently in 
progress. 

Assuming the local ordering pattern obtained in spin- 
wave theory, we can now investigate the question of long- 
range order by an Ising expansion for our original Hamil- 
tonian. We write the Hamiltonian as TC(X) = Hq + XH\, 
where, Hq consists of an Ising Hamiltonian, whose two 
ground states are ordered in the four-sublattice pattern 
and H\ = Ti — Hq, with TL being the full Hamiltonian 
in Eq. [j]. We develop expansions for the sublattice mag- 
netization in terms of the expansion parameter A. If the 
system has long-range order, the expansion should con- 
verge upto A = 1, whereas if it is disordered there should 
be a critical point at A < 1. Without explicit calcula- 
tions, we can make the following observation: because of 
the cancellation of the effective field from one chain to 
another, the purely one-dimensional graphs remain un- 
affected by the coupling between the chains. We know 
that these add up to zero sublattice magnetization at the 
Heisenberg point. The effects of the other graphs is to 
further reduce order, thus moving the Heisenberg system 
into a quantum disordered phase. We have carried out 
these expansions to order A 8 for a range of J'/ J values. 
The analysis suggests critical values less than unity, im- 
plying a disordered ground state for the Heisenberg limit . 

Additional insight can be gained by studying the J' 3> 
J limit of the model, where the Neel state on 45-degrees 
rotated lattice is stabilized by the J' interaction. As J 
increases from zero, quantum fluctuations grow stronger 
and at (J/ J') cr it — 0.76 (within linear spin-wave approx- 
imation) finally destroy long range order, analogous to 
Ji — Ji model j|) . That these two systems are very sim- 
ilar is also supported by the fact that the critical value 
of frustrating interaction J quoted above is exactly twice 
the corresponding critical value of J\ — J-i model. Given 
this analogy, it is tempting to speculate that disordered 
phase at J/ J' > 0.76 is a spontaneously dimerized one 
|l6| , ^2| , with vertical and horizontal chains formed by the 
exchange J being in the valence-bond-type state. How- 
ever, extensions of these results to J/ J' > 1 is prob- 
lematic as at J = J' this system classically has a finite 
ground state entropy and may change character at 
that point. Clearly more work is needed to understand 
the nature of the disordered phase. It is very interesting 
to ask what will happen upon introducing of a mobile 
charge carriers into the chains. We would like to point 
out apparent similarities of this problem with a very re- 
cent study of a square mesh of conducting horizontal and 
vertical stripes Cq] . 

In conclusion, we have identified a new frustrated 
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layered Heisenberg model, which could have a two- 
dimensional spin-liquid ground state. The magnetic 
properties of the materials Na2Ti2As20 and Na2Ti2Sb20 
maybe represented by such a model. 
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